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ABSTRACT 

We present a study of elemental abundances in a sample of thirteen Blue Compact Dwarf (BCD) 
galaxies, using the ~10-37/im high resolution spectra obtained with Spitzer/IRS. We derive the abun- 
dances of neon and sulfur for our sample using the infrared fine-structure lines probing regions which 
may be obscured by dust in the optical and compare our results with similar infrared studies of 
starburst galaxies from ISO. We find a good correlation between the neon and sulfur abundances, 
though sulfur is under-abundant relative to neon with respect to the solar value. A comparison of 
the elemental abundances (neon, sulfur) measured from the infrared data with those derived from the 
optical (neon, sulfur, oxygen) studies reveals a good overall agreement for sulfur, while the infrared 
derived neon abundances are slightly higher than the optical values. This indicates that either the 
metallicities of dust enshrouded regions in BCDs are similar to the optically accessible regions, or 
that if they are different they do not contribute substantially to the total infrared emission of the host 
galaxy 

Subject headings: dust, extinction — galaxies:starburst — abundances 



1. INTRODUCTION 

Blue Compact Dwarf Galaxies (BCDs) are dwarf 
galaxies with blue optical colors resulting from one or 
more intense bursts of star-formation, low luminosities 
{Mb > —18) and small sizes . The first BCD discovered 
was I Zw 18 bv lZwickvl ()1966[ ). which h ad the lowest oxy- 
gen a bundance observed in a galaxy (jSearle fc SargentI 
Il972| l. until the recent study of th e western component 
of SBS0335-052 (jlzotov et al.l[2n05l ). Although BCDs are 
defined mostly by their morphological parameters, they 
are globally found to have low heavy element abundances 
as measured from their HII regions (l/3OZ0'^l/2Zo). 
The low metallicity of BCDs is suggestive of a young age 
since their interstellar medium is chemically unevolved. 
However, some BCDs do display an older stellar popula- 
tion and have formed a large fraction of the i r stars more 
than IGyr ago (see iLoose fc ThuanI Il985t lAloisi et"aLl 
|2007( ). The plausible scenario that BCDs are young is in- 
triguing within the context of Cold Dark Matter models 
which predict that low-mass dwarf galaxies, originating 
from density perturbations much less massive than those 
producing the larger structures, can still be forming at 
the current epoch. However, despite the great success 
in detecting galaxies at high redshift over the past few 
years, bona fide young galaxies still remain extremely 
difficult to find in the local universe (jKunth fc Sargenti 
[TOSllKunth fc Ostlinl[2000l : lMadden et al.ll2006D . This is 
likely due to the observational bias of sampling mostly lu- 
minous more evolved galaxies at high redshifts. If some 
BCDs are truly young galaxies, they would provide an 
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ideal local laboratory to understand the galaxy forma- 
tion processes in the early universe. 

Over the past two decades, BCDs have been stud- 
ied extensively in many wavelengths using ground- 
based and space-born instruments (for a review see 
iKunth fc OsTlhH [200l . In the FUV, the Far Ultra- 
violet Spectroscopic Explorer (FUSE) has been used 
to study the che mical abundances i n the neutral gas 
in several BCDs (iThuan et all [200l lAloisi et al.] l2003t 
iLebouteiller et al.l l2004f )" Optical spectra have been 
obtained for a large number of BCDs and display 
strong narrow emission lines resulting from the in- 
tensive star-f ormation processes that take place in 
these systems (llzotov et al.lll997l:IIzotov fc Thuanlll999bl: 
iPustilnik et al.l 120051: iSalzer et all I2005D . The In- 
frared Space Observatory (ISO) revealed unexpect- 
edly that despite their low metallicities, BCDs, such 
as SBS 0335-052E, could still have c o pious emission 
from dust grains (IThuan et al.l fl999l: Madded I2OOOI : 
iMadden et all l2006t iPlante fc Sauvad I2002D. More 
recen tly, the Spitzer Space Telescope (iWerner et al.l 
l2004f ) has been used to observe these metal-poor 
dwarf systems in order to study their dust contin- 
uum properties a nd the polycyclic aro r natic hydrocarbon 
(PAH) features fHouck et al.' '2004b' 'Hogg e t all [20051: 
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Engelbracht ct al. 2005: Rosenberg ct al. 2006,; 
200ft lO'Halloran et al.ll200(THHunt et al.ll2006l : 
2007D . Finally, radio observations have also been per- 
formed for severa l BCDs to study t heir HI kinematics 
and distribution (|Thuan et al.l 1200^ and th ermal /non- 
therm al continuum emission properties (iHunt et al.l 
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Metallicity is a key parameter that influences the 
formation and evolution of both stars and galax- 
ies. Detailed studies of the elemental abundances of 
BCDs have already bee n carried out by seve r al groups 
(llzot ov fc ThuanI |i999b; Kniaze v et alT l2003t iShi et all 
120051) and the well known mctallicity-luminosity relation 
has also been studied in detail in the environment of 
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TABLE 1 

Observing Parameters of the Sample 



Object Name 


RA 


Doc 


AORKEY 


Observation 


Redshift 


On-source 


Time (Sen) 




(J2000) 


(J2000) 




Date 




SH 


LR 


TTarnI 1 






Qnn7i 04 


9004 07 1 7 


090fi 


480 


240 


IN v_i V_/ ± ±'±\J 




1 flrlni mArie 
-lUClUllIl^Ub 




9004 01 07 


OO^iO 


4S0 


240 


O J-J tJ L/ iJ iJ !j L/ tJ J-J 


n3h37m44 (Is 


_n^rin2Tn4ns 


1 1 7fi9S^fi 


2004-09-01 


0.0135 


1440 


960 






'JrtLJ. W 111 O T> O 


qnm 984 


2004-03-01 


~ 


480 


240 


TT7w4n 

i. J.ZJ W '-t\J 


^'^ll'^^m49 fit; 


v* O U. ^ O lllO ^ O 


Qnn7fii fi 


9004 0*^ 01 


009fi 


480 


240 




VJOli-LOlliXO . Uo 


ij u. u I? iiio I? o 


1 9n7fin'^9 


9004 1 9*^ 


0.0015 


120 


56 








1 9fi9fifiS;^ 


9004 11 11 




120 


56 


TZwl S 




t^OLJ. l_rtlll^OO 


qnr)8fi4n 


2004-03-27 


0025 


480 


240 








1 fi2n^^fi8 


2005-1 2-1 6 




2880 


1440 






i ou-oyiiioyo 


Qnn'^s9zi 

yuutjoz'i 


9004 1 9 OQ 





480 


240 


iVl i is. J. '±ij u 


1 1 Vi'^J^m'i'^ fit; 




9011712 


9004 19 19 

ZWUrt -LZ _LZ 


00*^9 


480 


240 




J_ lllu J-lilOO . Ob 


n9rl99TTi99Q 




900^1 01 O'i 




4S0 


240 








16204032 


2006-01-14 




1440 




SBS1210+537A 


12hl2m55.9s 


+53d27m38s 


8989952 


2004-06-06 




480 


240 


Toll214-277 


12hl7ml7.1s 


-28d02m33s 


9008128 


2004-06-28 


0.0260 


480 


240 


Tol65 


12h25m46.9s 


-36dl4m01s 


4829696 


2004-01-07 


0.0090 


480 


240 


UGCA292 


12h38m40.0s 


+32d46m01s 


4831232 


2004-01-07 


0.0010 


480 


240 


Toll304-353 


13h07m37.5s 


-35d38ml9s 


9006848 


2004-06-25 


0.0140 


480 


240 


Poxl86 


13h25m48.6s 


-Ild37m38s 


9007360 


2004-07-14 


0.0039 


480 


240 


CG0563 


14h52m05.7s 


+38dl0m59s 


8992512 


2005-05-30 


0.0324 


240 


120 


CG0598 


14h59m20.6s 


+42dl6ml0s 


8992256 


2005-03-19 


0.0575 


480 


240 


CG0752 


15h31m21.3s 


+47d01m24s 


8991744 


2005-03-19 


0.0211 


480 


240 


Mrkl499 


16h35m21.1s 


+52dl2m53s 


9011456 


2004-06-05 


0.0090 


480 


240 


[RC2]A2228-00 


22h30m33.9s 


-00d07m35s 


9006080 


2004-06-24 


0.0052 


480 


240 



Note. — The coordinates and rcdshifts of the objects are cited from the NASA/IPAC Extragalactic Database (NED). 
In this paper, we only include the analysis of thirteen out of twenty-two sources which have SNRs sufficient for our 
abundance study. CG0563, CG0598 an d CG0752 are included in the original sample as BCD candidates, however, they 
appear to be more starburst like (see iHao et al.ll2007f) . Thus even though they have high SNR, we do not include these 
three sources in this study. 



dwarf galaxies fSkillman ct al. 1 9891 : iHunter &: HofFmaiil 
\T999i .Melbourne fc Salzer. .20021 ). However, because 
these studies were performed in the optical, they were 
limited by the fact that the properties of some of the 
deeply obscured regions in the star- forming galaxies may 
remain inaccessibl e due to dust extinction at these wave- 
lengths. In fact, iThuan et all (jl999). using ISO, have 
shown that the eastern component of SBS 0335-052 does 
have an embedded super star cluster (SSC) that is in- 
visible in the optical while co ntributing ^^75% to the 
bolometric lu minosity (see also Plan t e fc Sauvag 
iHouck et all [2004b) . even though it has very low metal- 
licity (12-|-log(0/H)=7.33), which would in principle im- 
ply a low dust content. In addition to probing the dust 
enshrouded regions, emission in the infrared also has the 
advantage that the lines accessible at these wavelengths 
are less sensitive to the electron temperature fluctua- 
tions than the corresponding optical lines of the same 
ion. In the infrared, more ionization stages of an el- 
ement become available as well. The imp roved sensi- 
tivity of the Infrared Spectrograph (IRS^) f Houck et al.l 
[2004a) on Spitzer has enabled us to obtain for the first 
time infrared spectra for a much larger sample of BCDs 
than was previously possible ([Thu an et al."1999^. 'MaddenI 
120001 : IVerma et al.li2003l : iMartm- Hcrnandcz ct al. 2006), 
thus motivating this study to probe the heavy element 
abundances in BCDs. 

In this paper, we analyze Spitze/IRS spectra of thir- 
teen BCDs and present elemental abundances of neon 

^ The IRS was a collaborative venture between Cornell Univer- 
sity and Ball Aerospace Corporation funded by NASA through the 
Jet Propulsion Laboratory and the Ames Research Center. 



and sulfur, which are both primary elements produced 
by the same massive stars in the nuclear synthesis pro- 
cesses. In section 2, we describe the sample selection, ob- 
servations and data reduction. We present our results on 
the chemical abundances in section 3, followed by a com- 
parison of the optical and infrared derived abundances in 
section 4. We show the interplay between the abundances 
and PAH emission in section 5. Finally, we summarize 
our conclusions in section 6. 

2. OBSERVATIONS AND DATA REDUCTION 

As part of the IRS Guaranteed Time Observation 
(GTO) program, we have compiled a large sample of 
BCD candidates selected from the Secon d Byurakan Sur- 
vey (SBS), Bootes void galaxies (Kirs hner et al.l [19811 : 
[Popescu fc HodpI l2000f ). and other commonly studied 
BCDs. Details on the low-resolution spectra of the sam- 
ple have been published by Wu et al. (2006). 

We acquired the targets using the red (22/im) IRS 
peak-up camera in high accuracy mode to locate the mid- 
IR centroid of the source and then offset to the appropri- 
ate slit using the standard IRS staring observing mode. 
Subsequently, we obtained the 10— 37/im spectra for our 
sources using the ffiS' Short-High (SH, 9.9-19.6 /xm) and 
Long-high (LH, 18.7-37.2 /im) modules. No spectra of the 
background were obtained. The AORkey and on-source 
integration time of the twenty-two sources for which we 
obtain high-resolution spectroscopy are given in Table 
[TJ In this paper, we focus on the description and anal- 
ysis of thirteen BCDs from our program which have a 
signal-to-noise ratio (SNR) high enough for our elemen- 
tal abundance study. 

The data were processed by the Spitzer Science Cen- 
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ter (SSC) pipeline version 13.2. The three-dimensional 
data cubes were converted to two-dimensional slope im- 
ages after linearization correction, subtraction of darks, 
and cosmic-ray removal. The reduction of the spec- 
tral data started from intermediate pipeline products 
"droop" files, which only lacked stray light and flat- 
field correction. Individual pointings to each nod po- 
sition of the slit were co-added. The data from SH 
and LH were extracted using the full slit extraction 
method of a scri pt version of IR S data analysis pack- 
age SMART ( Hi gdon et all I2004D from the median of 
the combined images. The 1-dimensional spectra were 
flux calibrated by multiplying by a relative spectral re- 
sponse function (RSRF), which was created from the IRS 
standard star, £ Dra, for wh i ch accurate templates were 
available (jCohen et all 120031 : iSloan et al.|[200 7). As a fi- 
nal cosmetic step the ends of each order where the noise 
increases significantly were manually clipped. No scaling 
was needed between the adjacent orders within the same 
module. For faint sources, the spectra taken at nod po- 
sition 1 are severely affected by fringing problems, thus 
we only use the data at nod position 2 for these sources 
in our study. 

The fine-structure lines, [SIV] A10.51/im, 
[Nell] A12.8lAim, [Nelll] A15.55/imand [SIII] AAlS.Tl^m, 
33.48/im are clearly present in the majority of the BCDs 
in our sample^. We measure the line fluxes by fitting 
them with a Gaussian profile. Even though the [SIII] 
line is visible in both SH and LH we only use the 
18.71 /im line for deriving the ionic abundance. This is 
because by using lines within the same module (SH), 
we remove uncertainties due to the different sizes of 
the SH and LH slits. Furthermore, the [SIII] A33.48^m 
line is near the cut off of the LH module where the 
sensitivity drops dramatically, making its measurement 
more uncertain in some cases. 

3. ELEMENTAL ABUNDANCES OF BCDS 

In this section, we derive the neon and sulfur 
abundances using the new infrared data of our BCD 
sample. The elemental abu ndances of the compact 
HII r egions in our Gala xy (jMartin-Hernandez et al.l 
120021: iSimpson et al.l [200l . in the Magellanic Clouds 
(I Vermeil et al.l l2002f). as well as in starburst galax- 
ies (jVerma et al.l 120031: iMartin-Hernandez et all 120061) 
have already been studied extensivel y in the past usin g 
ISO, and more recently with Spitzer (jRubin et al.ll2007t ). 
There are several metho ds to derive the chemical abun- 
dances (for a review see IStasihska 1 120071) . Here we use 
an empirical method, which derives ionic abundances di- 
rectly from the observed lines of the relavent ions. To 
do so, we need to have the flux of at least one hydro- 
gen recombination line, the dust extinction, as well as 
the electron temperature and density of the interstellar 
medium (ISM). 



3.1. Electron density and temperature 

The electron density (Ne) could in principle be 
determined by comparing the measured ratio of 
[SIII]18.71^m/33.48^m to the expect ed theoretical valu e 
using the corresponding 5'-curve (see iHouck et al]|1984D . 

6 See Fig. 3 of IWu et all l l200l) . or IHao et all | |2007|) for the 
reduced spectra. 



However, the values for the ratio are in the horizontal 
part of the 5'-curve and as a result we cannot use these 
two lines to accurately constrain the electron density of 
these systems. Furthermore, the 33.48/im [SIII] line is 
located at the edge of the LH slit, where the sensitivity 
drops dramatically, and considering that these emission 
lines are weak, the measured line flux for the 33.48 /xm 
[SIII] has a large uncertainty. Since the infrared determi- 
nation of elemental abundances does not depend strongly 
on the density, we adopted the optically derived elec- 
tron densities from the literature, which range from 10 
to 3000 cm~'^ (see Tabled]). For sources where such in- 
formation is not available, we adopt a typical electron 
density of 100 cm~'^ in this paper. 

Since BCDs in general have low metallicities, it is 
expected that BCDs would have relatively high elec- 
tron temperatures (Tg). Fortunately, the infrared 
lines are much less sensitive to the uncertainties in 
the electron temperature compared to the optical 
(jBernard Salas et al.ll200ll ). As a result we adopt elec- 
tron temperatures derived from the optical studies found 
in the literature. For sources which do not have di- 
rect measurements (which happen to be high metallicity 
sources in our sample, e.g. NGC1140, UGC4274), we use 
a representative temperature of Tg ~ 10,000 K (see Table 

3.2. Ionized hydrogen flux estimates 

In order to derive the ionic abundances, we also need 
to obtain information on at least one hydrogen recom- 
bination line (usually H/3). For sources in which the 
Humphreys a (Huq, 12.37 /im) line is detected, we use 
it to convert to H/3. From the extinction corrected Hua 
flux, we esti mate the H/3 flux us i ng the values given in 
the tables of iHummer fc Storevl (119871 ) for Case B re- 
combination. This has the advantage that we do not 
need to correct for aperture effects because the neon and 
sulfur lines we use to derive the abundances also reside 
in the same module (SH). The 12.37/Ltm line is gener- 
ally very weak and it is detected in six of our thirteen 
sources. For SBS0335-052E and IZwl8, which have ra- 
dio continuum data, we calculate the H/? flux from the 
thermal free-free emission using the fraction included in 
the aperture of our SH slit at the time of the observa- 
tion, and no correction for extinction is needed. For the 
rema ining four source s which have Ha images available 
from lGil de Paz et al.l (|2003f ). we overlay the SH slit on 
each Ha image to calculate what fraction of the Ha emis- 
sion is included in our slit. Then we derive the corre- 
sponding H/3 flux inside the SH slit from the extinction- 
corrected Ha flux. Finally, for UM461, we use the H/3 
flux of the whole galaxy from the integrated spectra 
(jMoustakas fc Kennicut?l2006f) . UM461 is a very com- 
pact source, thus fully encompassed within the SH sht. 

3.3. Extinction correction 

All our measurements, both the fluxes from the Ha im- 
ages, as well as the infrared lines, have been corrected for 
extinction. From the low-resolution spectra of our BCDs 
(IWu et al.ll200l . we find that most of the sources do not 
show a strong 9.7 /im silicate feature, thus indicating an 
intrinsically relatively low dust extinction. Due to the 
faintncss of our targets at A <9.7/im, the determination 
of the continuum on the blue side of the 9.7 /im feature is 
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TABLE 2 
Optical Properties of the Sources 



Object F(H/3) (xlO'^'^ergs s'^cm'^) Eg^v Te Ne ref 

Hua-derived Ha-derived radio-derived optical mag K cm~'^ 



Haroll 


252±40 




0.41 


13700 


10 


NGC1140 


67±9 




0.10 


10000 


100 


SBS0335-052E 


<22.4 




6.5±0.06 ■ ■ ■ 2.76 


20000 


200 


NGC1569 


393±19 




0.65 


12000 


100 


IIZw40 


363±10 


282±22 


0.79 


13000 


190 


UGC4274 


<59.6 


11.5±1.0 





10000 


62 


IZwl8 


<12.6 


8.0±0.6 


6.1±0.6 •■■ 0.08 


19000 


100 


VIIZw403 


<24.6 


11.8±1.2 





14800 


100 


Mrkl450 


<14.8 


15.8±1.4 


0.10 


12500 


100 


UM461 


<21.0 




13.6±0.6 0.08 


16100 


200 


Toll214-277 


23.3±2.8 




0.03 


19790 


400 


Tol65 


<18.2 


9.3±0.4 


0.08 


17320 


50 


Mrkl499 


14.6±1.9 


15.1±1.3 


0.17 


12600 


3267 



(1) 

(2)- 
(3),(4),{5) 
(6) 
(2),{7) 
{7),(8)- 
{7),(9),{10) 
(7),{11) 
(7),{12) 
(13),(14) 

(15) 
(7),{15) 
(7),{16) 



Ref erenc es. — (1) BorKvall & Oumlstlin (2002), (2) Guscva ot al. (2000), (3) Izotov ct al. (2006), (4) Hunt c t al. 
1 2004[) . (51 IHou ck et al. |2004b), (6) Kobulnickv & Skillman (1997), (7) Gil dc Paz ct al. (2003), (8) Ho ct al. 
a997l). (9) ICann on ct al. (2005), (10) Izotov ct al. (1999a), (11) Izotov ct al. (1997), (12) Izotov ct al. (1994|). (13) 
IMoustakas fc Kcn nicutt (2006), ( 14) Izotov &: Thuag (,1998.). (15) .Izotov ct al. (2001) , (16) Kong ct al.. (,2002f) 

Note, — When no data are available in the literature, the electron temperature and densities are assumed to be 
10 000 K ( for NGC1140 and UGC4274) and 100 cm"'' (for NGC1140 and NGC1569) respectively. When the Hu a line 
is detected, the H/3 flux is preferentially derived from this line. The upper limits are listed for the non-detections. For 
the remaining sources, wc derived the ¥Lj3 flux either from the thermal component of the radio continuum, or from the 
extinction corrected Ha flux inside the SH slit. For UM461, the H/3 flux is for the whole galaxy from the integrated 
optical spectra. 



often poor and the uncertainties in estimating the mid- 
IR extinction from this are rather large (see discussion in 
ISpoon et al.l[200l . Thus for sources where optical spec- 
tra are available, we adopt the Eb-v values calculated 
from the hydrogen recombination lines in the literature. 
Only for SBS 0335-052E, which has a high quality mid-IR 
spectrum and evidence for an embedded SSC, the opti- 
cal Es-v is not a good estimate of its extinction, thus 
we use the exti nction estimated fro m the depth of the 
silicate feature (jHouck e t alJ l2 004bl). T hroughout this 
paper, we adopt the [M uks et all (|1994l ) extinction law, 
though using the iDraing (j2003) law would produce very 
similar results. The Eb-v magnitudes are provided in 
Table [21 

3.4. Neon and sulfur abundances determination 

We use the electron density and temperature given 
in Table [2l together with the fine structure lines de- 
tected in the IRS high-resolution spectra, to calcu- 
late the abundances of neon and sulfur. The equation 
used to determine the ion ic abundance is described in 
[Bernard Salas et al] (|200lD and it is: 

Nion _ lion T. J Aul QHg / , . 

Np " /h„ 'Ah, ^ui \NionJ ^ ' 

where Iion/lH/3 is the measured flux of the fine- 
structure line normalized to H/3; Np is the density of 
the ionized hydrogen; Aui and Ah/3 are the wavelengths 
of the line and H/3; an/s is the effective recombination 
coefficient for H/3; Aui is the Einstein spontaneous tran- 
sition rate for the line and Nu/Nion is the ratio of the 
population of the level from which the line originates 
to the total population of the ion. This ratio is deter- 
mined by solving the statistical equilibrium equation for 
a five level system and norm alizing the total number of 
ions to be unity l|Osterbrock ,1989 ). The effective col- 
lisional strengths used to derive the population of lev- 



els were obtai ned from the appropr iate reference of the 
IRON project ([Hummer et al.lll993f )'^. 

The most important ionization stages of neon and 
sulfur are available in the infrared. We detected 
[SIV] A10.51 /im and [Nelll] A15.56/L(m in aU of our BCDs 
while the [NeH] A12.81 fim line is detected in nine ob- 
jects and [SIII] A18.71 /xm in twelve sources (see Table 
[21). For su lfur, SH has i t s stro ngest emission lines in 
the optical. [Vermeil et al.l ([2002[ ) have shown that SH is 
typically less than 10% of SIII in the Magellanic cloud 
HII regions they studied. Because BCDs are typically 
high-excitation objects, we do not expect to have a sig- 
nificant contribution from SII. The optical studies reveal 
that the ioni c abundance of SII is typically 10—20% of 
that of SIII ([Izotov et al.[[1994 [l997.) , thus we keep in 
mind that this might result in a ~10% underestimate in 
our infrared derived sulfur abundance. This is addressed 
further in the following subsection where we discuss the 
properties of the individual objects. The presence of the 
[OIV] A25.89/im line in most of our sources (except for 
VIIZw403, UM461, Tol65 and Mrkl499) indicates that 
some [NelV] might be present, even though the [OIV] 
line is typically weak. For example, in Toll214-277, 
where [OIV] is detected, the fiux of [NeIV]A4725 from 
optical spectra is less than 1% of that of [NeIII]A3868 
([Izotov et al.|[200"^ . indicat ing that the contribu t ion o f 
[NelV] is not significant. [Tsamis fc Peq uignol^ (| 2005[ ) 
also have shown in their photoionization modelling study 
of 30 Doradus that the contribution from NelV to the 
total neon abundance is much less than 1%. As a re- 
sult, rather than applying ionization correction factors 
(ICFs) which could introduce an unknown systematic un- 
certainty in our study, we sum the corresponding ionic 
abundances for the most important ionization stages to 
determine the total element abundance and the results 

^ http:/ /www.usm.uni-muenchen.de/people/ip/iron- 

project.html 
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TABLE 3 

Fine-structure Line Fluxes and Ionic Abundances 



Line Fluxes ( xlO -'^''ergs cm ^ s ^) 



Ionic Abundances (xlO ^) 



Object 


SI V 10 Slum 


[Ncil](l2.8l/^m) 


[Nelll](l5.55/im) 


Sill 18 71itm 




Nm tt / 

^^NcII/ ^^p 


^^NcIII/ ^^p 


Ncttt/N^ 
^^SIII/ ^^P 


Haroll 


43.95±0.14 


31.23±0.23 


98.14±0.65 


45.24±0.25 


0.353±0.056 


15.4±2.4 


23.2±3.7 


1.89±0.30 


NGC1140 


12.14±0.13 


11.41±0.20 


37.79±0.32 


20.04dz0.11 


0.375±0.051 


22.5±3.0 


35.5±4.8 


3.56±0.48 


SBS0335-052E 


1.56±0.03 


<0.14 


1.43±0.02 


0.43±0.03 


0.527±0.011 


<2.27 


11.2±1.6 


0.535±0.038 


NGC1569 


147.26±0.55 


15.70±0.17 


175.85±0.59 


72.70±0.23 


0.817±0.040 


5.28±0.26 


28.2±1.4 


2.12±0.10 


IIZw40 


185.72±1.73 


6.24±0.11 


112.65±0.83 


45.80±0.21 


1.12±0.03 


2.23±0.07 


19.3±0.6 


1.37±0.04 


UGC4274 


4.18±0.12 


8.87±0.11 


12.84±0.12 


12.63±0.11 


0.743±0.068 


102dz9 


70.3±6.1 


13.2±1.2 


IZwlS 


0.48±0.03 


0.09±0.01 


0.46±0.02 


0.23±0.02 


0.125±0.015 


1.41±0.21 


3.56±0.38 


0.276±0.036 


VIIZw403 


0.99±0.05 


<0.142 


0.88±0.02 


0.87±0.02 


0.311±0.035 


<2.79 


8.41±0.88 


1.39±0.14 


Mrkl450 


7.34±0.04 


1.34±0.03 


9.61±0.04 


4.37±0.03 


0.928±0.223 


10.70±2.57 


36.9±8.9 


2.98±0.71 


UM461 


4.58±0.04 


0.16±0.02 


2.83±0.03 


0.812±0.03 


0.636±0.029 


1.38±0.18 


11.8±0.5 


0.563±0.032 


Toll214-277 


0.88±0.02 


<0.10 


0.61±0.02 


<0.19 


0.061±0.007 


<0.413 


1.23±0.15 


<0.059 


Tol65 


0.69±0.02 


<0.14 


0.90±0.02 


0.32±0.02 


0.129±0.007 


<1.62 


5.07±0.25 


0.293±0.022 


Mrkl499 


1.74±0.02 


1.88±0.03 


5.03±0.03 


2.91±0.02 


0.286±0.037 


16.5±2.2 


21.6±2.8 


2.28±0.30 


Note. — The observed line fluxes 


are measured from the IRS high resolution spectra of the sources. Background 


emission has not been subtracted, but this 



does not affect the flux of the flne-structure lines. 




Ne"/Ne* 

Fig. 1.— A plot of Ne++/Ne+ vs S^+/S++, which traces the 
hardness of the radiation field. The solid line is a fit to the data. 



are shown in Table [3l 

Two successive stages of ionization of a given element 
(X) can be used to measure the state of the ionization of 
the ISM, which, on first order, depends on the ionization 
parameter , U, and the hardness of the ionizing radia- 
tion ( Vilc hez &: Pag el 1988). For a given U, the ratio 
X+*+^/X+* is an indicator of the hardness of the stellar 
radiation field. Here we use the ionic pairs of Ne"'"+/Ne+ 
and S"'"'^/S"'"+ as indicators of the hardness of the ion- 
ization field and plot them in Fig. 1. It is clear from 
the figure that these two ratios increase proportionally. 
This indicates that the hardening of the radiation field 
affects similarly the various elements in the full range of 
the ionizing continua. 

The major uncertainty in our study of the elemental 
abundances originates from the II/3 fluxes used as refer- 
ences. For sources where the 12.37/im Hua line is de- 
tected, we have no uncertainties from matching the slit 
apertures, and the extinction correction is very small. 
For sources where we estimate the EE/? flux from the ther- 
mal radio continuum emission, we are free of extinction 
corrections, even though we still need to properly account 
for the radio emission which falls within the IRS slit. For 
BCDs with no Hua emission or radio data available in 



the literature, we estimate the H/3 flux using archival 
Ha images of the objects. This is clearly more challeng- 
ing because in addition to the aperture corrections, the 
extinction correction is higher which could introduce a 
systematic uncertainty as large as 25%. In the following 
analysis, we indicate the sources with or without the de- 
tection of Hu a line using different symbols on the plots 
(see Fig. 2, 3 and 4). The measurement uncertainties in 
the infrared lines are typically only a few percent, but 
could grow to ^^10% in the case of fainter sources. The 
uncertainties we list in Table [3] only account for measure- 
ment uncertainties. 

3.5. Individual Objects 

Some details on the individual objects and how we 
measured the necessary parameters are provided in the 
following paragraphs. 

Haroll: The 12.37/im Hua line is clearly detected 
in the SH spectrum. The source has multiple nucle i, 
designated as A, B and C (|Bergvall fc Oumlsthn|[200l) . 
The optical spectra from which the oxygen abundance 
was derived corresponds to the central region. Because 
all the lines we are using are from the SH spectrum, we 
are relatively free from extinction problems and there is 
no need for aperture correction. The contribution of SII 
may add ^10% to the total sulfur abundance. 

NGC1140: As with Haroll, the 12.37/im Hua line 
is also clearly detected in the SH spectrum of NGC1140. 
No measurement was found in the literature for the elec- 
tron temperature or density for this source, thus we as- 
sume 10000 K and 100 cm^'^ respectively, which are typi- 
cal for BCDs. However, because its oxygen abundance is 
more than half solar, the electron temperature could also 
be lower. If we use a Tg of 5000 K and re-derive the abun- 
dances, we find that the neon abundance would increase 
from 5.8x 10^^ to 8.5xl0~^ and the sulfur abundance 
from 3.9x10-6 to 5.5 xlO'^. The contribution of SII may 
add ^^10% to the total sulfur abundance. 

SBS0335-052E: We do not detect the 12.37/im Hua 
line in the SH spectrum. We use th e thermal compon ent 
of the 5 GHz radio continuum from lHunt et al.l ()2004f ) to 
convert to the H/3 flux, and thus no extinction correc- 
tion is needed. [Nell] is not detected and we provide the 
upper limit in Table [3l The ionic abundance of SII is 
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TABLE 4 
Elemental Abundances 



Object Nc/H(xlO-'^) S/H (xlQ-*^) O/U (xlQ--') Nc (Zq) S (Zq) O (Zq) ref (O/H) 



Haroll 


38.6±4.4 


2.24±0.31 


7.94 


0.33±0.04 


0.16±0.02 


0.16 


(1) 


NGC1140 


58.0±5.7 


3.94±0.48 


28.8 


0.48±0.05 


0.28±0.03 


0.59 


(2) 


SBS0335-052E 


[11.2, 13.5]=' 


1.06±0.04 


1.95±0.09 


[0.09, 0.11]=" 


0.08±0.003 


0.03±0.001 


(3) 


NGC1569 


33.5±1.4 


2.94±0.11 


15.8±0.5 


0.28±0.01 


0.21±0.01 


0.32±0.01 


(4) 


IIZw40 


21.5±0.6 


2.49±0.05 


12.2±0.4 


0.17±0.005 


0.18±0.004 


0.25±0.01 


(2) 


UGC4274 


172.3±10.9 


13.94±1.20 


33.3 


1.4±0.1 


l.OiO.l 


0.68 


(5) 


IZwlS 


5.0±0.4 


0.40±0.04 


1.47±0.09 


0.04±0.003 


0.03±0.003 


0.03±0.002 


(6) 


VIIZw403 


[8.4, 11. 2]=' 


1.70±0.14 


4.90±0.11 


[0.07, 0.09]=» 


0.12±0.01 


0.10±0.002 


(7) 


Mrkl450 


47.6±9.2 


3.98±0.74 


9.55±0.11 


0.40±0.08 


0.28±0.06 


0.19±0.002 


(8) 


UM461 


13.2±0.5 


1.20±0.04 


6.10±0.40 


O.lliO.Ol 


0.09±0.003 


0.12±0.01 


(9) 


Toll214-277 


[1.2, 1.6f 


[0.06,0.12]'' 


3.45±0.10 


[0.010, 0.013]=» 


[0.004, 0.008]*= 


0.07±0.002 


(10) 


Tol65 


[5.1, 6.7]=' 


0.42±0.07 


3.48±0.10 


[0.04,0.06]=^ 


0.03±0.002 


0.07±0.02 


(10) 


Mrkl499 


38.1±3.6 


2.57±0.30 


13.2 


0.32±0.03 


0.18±0.02 


0.27 


(11) 



Ref eren ces. — ( 1) LB crgvall fc Oumlstlin I 12003 ). (2 ) IGuseva et al.l ||2000D . (3) llzotov et al.l ||2006D. (4) [ Kobulnicky fe SkillmanI 
Ill997f). fSllHo et al.l 11997). (6) Izotov ct al ll999aD . (7) IIzotov et aLllll997l) . (8) IIzotov et al.llll994D . (9) ITzotov &: ThuanI ||1998|) . (10) 
llzotov et al.l II2001I) . (11) iKong ct al, (i200a) 

Note. — We adopt the following values for the solar abundances: (Nc/H)0 = 1.2 X lO"**, (S/H)0 = 1.4x lO"'^ and (O/H)0=4.6x 10"* 
(See section 4.1). 

^ Wc provide a range for the neon abundance by taking as the lower value the ionic abundance of Nelll, while to obtain the upper 
value wc add the upper limit for Nell.^ We provide a range for the sulfur by taking as the lower value the ionic abundance of SIV, 
while to obtain the upper value wc add the upper limit of SIII. 



~25% of the SIII determined from the optical study of 
llzotov et alj ()2006f ) and may add 12% to the total sulfur 
abundance. We also find higher neon and sulfur abun- 
dances compared wit h the oxygen and possible implica- 
tions are discussed in lHouck et al.l (j2004b). 

NGC1569: The 12.37 /im Una line is clearly detected 
in the SH spectrum, thus providing a direct estimate of 
the H/3 flux inside the SH slit. This galaxy is extended 
in the mid-IR (see Fig. 4 in IWu et al.l I2006D and op- 
tical spe ctra have been taken for seve ral of the bright 
knots bv lKobulnickv fc SkillmanI (|199 7^ in order to study 
the chemical gradient and inhomogeneities of the source. 
These authors found very little variation in the metallic- 
ity. Our infrared-derived abundances are in rough agree- 
ment with the optical results. The ionic abundance of 
SII is ~13% that o f the SIII estimated from the optical 
study (|Kobulnickv fc Skillman 1997 ) and may add 9% to 
the total sulfur abundance. 

IIZw40: We use the 12.37/im Hua line to convert to 
up flux. The derived neon and sulfur abundances agree 
with each other but appear to be lower than oxygen with 
respect to the solar v alues. If we estimate t he H/3 flux 
from the Ha image of lGil de Paz~iF^ (p003h . we find a 
value ^22% lower than the first, which would increase the 
neon and sulfur abundances by 22% accordingly. This 
suggests, that when using the H/3 flux derived from Ha 
image, the systematic error in our measurement is prob- 
ably no better than 25%. The ionic abun dance of SII is 
~13% that of SIII from the optical study of lGuseva et al.l 
dlOOOf) and may add 7% to the total sulfur abundance. 

UGC 4274: We use an H/3 fl ux derived from the Ha 
image of lGiTde Paz et all (pOOSi ). Similarly to NGC1569, 
the galaxy is extended and the peak of the infrared cen- 
troid is displaced from the optical peak position. The 
derived neon and sulfur abundances are both super solar. 
The oxygen abundance is not directly available from the 
literature, thus we use the [Nil] /Ha method ()Ho et al.l 
119971: iDenicolo et"aLl |2002|) to derive O/H, which bears 
a large uncertainty. The contribution of SII may add 



^10% to the total sulfur abundance. 

IZwlS: The H/3 flux is derived from the thermal com- 
ponent of the radio continuum inside the SH slit. The 
ionic abunda nce of SII is ^21% th at of SIII from the op- 
tical study of llzotov et al.l ()1999al ) and may add 15% to 
the total sulfur abundance. A more detailed d iscussion 
on this object can be found in lWu et"an ()2007| ). 

VIIZw4 03: The H/3 flux is der ived by using the Ha 
image from lGil de Paz et al.l ()2003f ). The optical and in- 
frared centroids do not overlap. The SH spectrum for 
this source is noisy and [Nell] A12.81 /im is not detected, 
thus we have only a lower limit on the neon abundance. 
The upper limit of the [Nell] line indicates that it could 
add less than 25% to the total elemental abundance of 
neon which is presented in Table |4l The ionic abun- 
d ance of SII is ^24% that of SIII from the optical work 
of llzotov et all ()1997D and may add 20% to the total sul- 
fur abundance. 

Mrkl450: As with VIIZw403, we use the hydrogen 
flux as given by the Ha flux included in the SH slit. 
Neon and sulfur abundances appear to be higher than 
the oxygen abundance with respect to the solar values. 
The ionic abun dance of SII is ~18% that of SIII from the 
optical study of llzotov et al.l (|1994D and may add 12% to 
the total sulfur abundance. 

UM461: The H/3 flux is estimated fr o m the inte- 
grated spectra of iMoustakas fc KennicuttI (j2006D , and 
the source is small enough to be fully included by the 
SH slit. The derived neon and sulfur abundances are in 
good agreement with oxygen metallicity from the optical. 
The ioni c abundance of SII is ^ 17% that of SIII from the 
optical (|Izotov fc ThuanI 119981 ) and may add 8% to the 
total sulfur abundance. 

Toll214-277: The 12.37 /um Hua line is marginally 
detected in the SH spectrum, and we use this line to 
derive the H/3 flux. The spectrum is very noisy and no 
[Nell] A12.81 ^m or [SIII] A18.71 /im lines are visible. We 
use the ionic abundances of [Nelll] and [SIV] as the lower 
limits of the neon and sulfur abundances. The upper 
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limit of [Nell] amounts to 33% of the ionic abundance of 
[Nelll] while the upper limit of [SIII] is nearly as much 
as the abundance of [SIV]. We give the range of their 
elemental abundances in Table IH The ionic abundance 
of SII is ^18% th at of SIII from the optical study of 
llzotov et alJ (|2001[ ) and may add ^9% to the total sulfur 
abundance. 

Tol65: We use the Ha image (|Gil de Paz et al.ll2003D 
to calculate the H/3 flux. [Nell] A12.81 fj,m is not detected 
in the SH spectrum, thus our derived neon abundance is a 
lower limit. The upper limit of Nell is ~32% of the ionic 
abundance of Nelll, and the range of neon abundance is 
provided in Table H) The metallicity of sulfur relative to 
solar is lower than oxygen. The ioni c abundance of SII 
is ~26% that of SIII from the optical ijlzotov et al.ll2001h 
and may add 19% to the total sulfur abundance. 

Mrkl499: The 12.37 ^m Hua line is detected in the 
SH spectrum. There was a bad pixel at the peak of this 
line and we interpolated its value using adjacent pixel val- 
ues. The infrared centroid appears to be slightly shifted 
from the optical one. The H/5 flux inside the SH slit de- 
rived from the Ha image is very similar to that derived 
from the Hua line (within ~10%). The contribution of 
SII may add ~10% to the total sulfur abundance. 

4. DISCUSSION 

In this section, we compare the abundances derived 
from the IR spectra with the results from the optical 
studies. 

4.1. Solar abundances 

The solar photospheric abundances have changed 
drastically over the past few years (see discussion in 
iPottasch fc Bernard-Salasll2006f ). In some cases, the so- 
lar abundance estimates for the same element from dif- 
ferent authors could differ by a factor of two. In this 
paper, when we quote the elemental abundance in so- 
lar units (e.g. Table [4]), we adopt the following values: 
(Ne/H)0=l. 2x10-4 from iFeldman fc WidineT /200a) . 

(l2005l) and 



Asplund et al.l 



Allende Prieto et al 




(S/H)0 = 1.4xlO-^ from 
(O/H)0=4.9xlO-4 from 

We should note though that Asplun d et al.l 
have also reported th at (Ne/H)Q=6.9x 10~^ and 
IGrevesse fc Sauval lHOOSh h ave given (S/H)0=2.1xlO-^ 
while lAnders fc Grevessd (| 19890 have reported that 
(O/H)0=8.5x 10~4. These values represent the extremes 
for the sol ar neon, sulfur and oxygeii abun dance deter- 
minations. IPottasch fc Bernard-Sala^ (|200 6^ have shown 
that the higher neon solar value is favored in their sam- 
ple of planetary nebulae (PNe) . In the present paper we 
indicate the range of solar abundances in our plots and 
we also discuss the impact of different solar values on our 
results. 

4.2. Neon and sulfur abundances 

As discussed in iThuan et al.l (|1995( ). it is assumed — 
and has been shown in the optical — that the a-element- 
to-oxygen abundance ratios do not vary with oxygen 
abundance. This is due to the fact that those elements 
are produced by the same massive stars (M>1OM0) re- 
sponsible for oxygen production. We test this assump- 
tion in the infrared by studying the abundances of neon 
and sulfur as derived from our IRS data. 




10" 



10" 



10" 



Fig. 2.— A plot of the Ne/H vs S/H abundance. Our BCDs are 
indicated by the filled circles if Hu a is detected and by diamonds 
if it is not detected. The sources that are shown by only the error 
bars are those that have non-detection of [Nell] or [SIII ] (See note 
of Table |4]l . The starburst galaxies from IVerma et al.l I I2003I1 are 
marked with the stars. We use the grey band to indicate the locus 
on the plot where the ratio of Ne/S would be consistent with the 
different solar values for neon and sulfur abundances. 




Ne/H 

Fig. 3. — The abundance ratio of Ne/S as a function of Ne/H. 
The symbols are the same as in Fig. 1. The grey area indicates 
the range for the values of the Ne/S ratio encountered in the solar 
neighborhood. 

In Fig. 2, we plot the abundances of ne on and sulfur 
as der ived from our high-resolution data. IVerma et al.l 
()2003l ) have found a positive correlation between the 
neon and argon abundances for their sample of starburst 
galaxies while their data show no correlation between the 
sulfur and neon and/or argon abundances (indicated as 
stars on Fig. 2). However, our sources show that the 
neon and sulfur abundances scale with each other. In 
the same figure, we also plot the proportionality line of 
the ratio for (Ne/S)0. The maximum and minimum val- 
ues of the solar neon and sulfur abundances are indicated 
by the width of the grey band, and we find that most of 
our BCDs have ratios above those values. 

In Fig. 3, we examine the ratio of Ne/S as a function of 
Ne/H derived from the IRS data. We see no correlation 
between the ratio of Ne/ S with respect to Ne/ H, in agree- 
ment with the results of iThuan et al.l ()1995t ). The aver- 
age Ne/S ratio we found for these BCDs is ^^11. 4. This 
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is similar with the Ne /S ratio of 14.3 in the Orion Nebula 
(|Simpson et al.ll2004f ). The Ne/S ratios found in the HII 
regions of M83 are higher (~24— 42) and d ecrease with 
incre asing deprojected galactocentric radii (jRubin et al.l 
|2007( ). These authors have also found that t he Ne/S ra- 
tios in M33 are -^12-21 (Rubin et al.ll2006f l. We notice 
that all the BCDs, as well as the HII regions in M33, M83 
and the Orion Nebula, have larger Ne/S ratios than those 
found in the solar neighborhood, while their Ne/S ratios 
from optical measurements are more consistent with the 
solar ratios. The higher Ne/S ratios from infrared stud- 
ies have already been found in other type of objects, in- 
cluding PNe, starburst galaxies, and other HII regions 
jMarigo et al. 2003; Vcrma ct al. 20031; iHenrv et all 
[2004t IPottasch fc Bernard-Salasi 12006 1. Some of these 
works point out that the solar sulfur abundance 
could have been underestimated and they also ap- 
pear to favor the higher values for the solar neon 
abundance. It could also be due to a differen- 
tial depletion of sulfur onto dust grai ns compared to 
neon as suggested by several studies (.Simpson fc Rubinl 
19901: IVerma et al.ll2003HPottasch fc Bernard-Salasl2006t 
Bernard Salas et al.ll2007f K 

4.3. Comparison with optically derived abundances 

As has been discussed earlier in this paper, one of 
the major advantages of using IR lines for estimating 
abundances is that one can probe emission from dust 
enshrouded regions without the uncertainties introduced 
from extinction corrections. In this subsection, we com- 
pare the newly derived neon and sulfur abundances of 
our BCD sample using the infrared data, with the abun- 
dances of the same elements as well as those for oxygen 
estimated from optical studies available in the literature. 

In order to determine elemental abundances, the op- 
tical studies usually adopt a two-zone photoionized HII 
region model: a high-ionization zone with temperature 
Te(OIII) and a low- ionization zone with temperature 
Te(OII). Because the infrared lines are much less sen- 
sitive to variations in Te (which affect the calculation 
on the population of the various atomic levels), in the 
subsequent calculations, we assume a constant electron 
temperature for our analysis. 

In Fig. 4 we compare the infrared derived neon and 
sulfur abundances of the BCD sample with their optical 
metallicities, using the oxygen abundances from the lit- 
erature (see also Table |4]). The infrared neon abundances 
with respect to solar are slightly higher than the oxygen 
(see Fig. 4a), but overall, there is a good agreement be- 
tween the infrared measurements and the optical results. 
This indicates that there is little dust enshrouded gas, or 
that these regions have similar metallicities, or that they 
do not contribute much to the total integrated emission 
of the galaxies. In SBS0335-052E, as well as Mrkl450, 
UGC4274 and Haroll, the infrared derived metallicity of 
neon is more than twice that of oxygen compared to the 
solar values. This could be due to the different regions we 
are probing using the infrared, which could be obscured 
in the optical. SBS0335-052E is known to have an em- 
bed ded super star cluste r that is invisible in the optical 
(see iHouck e"raIll2004bD . Both Haroll and UGC4274 
are somewhat extended, and detailed spectroscopy over 
the whole galaxy would be needed to directly compare 
the infrared results with the optical ones. For Mrkl450, 
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Fig. 4. — a)The Ne/H vs O/H abundances. The symbols are 
the same as in Fig. 1. The grey band represents the locus on 
the plot where sources of different metallicities would have a Ne/O 
intrinsic ratio similar to the values found in the solar neighborhood, 
b) Same as a), but for S/H vs O/H abundances. 

which is a point source for Spitzer/IRS, the abundance 
we derived using the infrared data is higher than the 
optical value and it is likely due to the fact that the op- 
tical estimates are more susceptible to the uncertainties 
in the electron temperature. If we lower Te from 12 500 K 
to 10 000 K, the oxygen abundance calculated from the 
optical data would double, while the neon abundance 
measured from the infrared would only increase by 5%. 
In Fig. 4b we show that the infrared derived sulfur abun- 
dances agree well with the oxygen abundance relative to 
solar. 

For nine out of the thirteen objects we studied, opti- 
cally derived neon and sulfur abundances are also avail- 
able in the literature. The optical data have matched 
apertures, in most cases an accurate Te measurement, 
a direct measure of the H/? flux, and reddening derived 
from the same data. Our infrared results are less sen- 
sitive to the electron temperature, but affected by the 
uncertainty in the H/3 flux if the Hua line is not de- 
tected and the uncertainty could be as large as 25%. 
In Fig. 5a, we plot the optically derived Ne/H against 
the infrared Ne/H. The solid line is the 1:1 proportion 
line for the infrared and optical derived abundances and 
we find that most of the sources are located above this 
line. A possible explanation for the higher infrared de- 
rived neon abundances could be the presence of dust 
enshrouded regions which might have higher heavy ele- 
ment abundances. This could explain the higher infrared 
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TABLE 5 

PAH MEASUREMENTS OF THE SAMPLE 



Object Name 




PAH EW (tim) 




Integrated Flux (xlO 


— 14 

ergs cm 


^ S ^ ) 




6.2 fim 


7.7 ^im 


8.6 fim 


11.2 /.tm 


6.2 fini 


7.7 fim 


8.6 fj,m 


11.2 fim 


Haroll 


0.120±0.003 


0.227±0.005 


0.042±0.003 


0.099±0.002 


78.3±1.8 


180±4 


25.4±1.6 


67.2±1.3 


NGC1140 


0.479±0.033 


0.568±0.024 


0.091±0.006 


0.525±0.032 


42.1±1.4 


78.3±1.8 


11.1±0.7 


50.2±2.5 


SBS0335-052E 


<0.035 






<0.015 


<1.7 






<1.0 


NGC1569 


0.232±0.013 


0.380±0.019 


0.007±0.005 


0.129±0.008 


50.4±2.4 


118±6 


2.7±1.6 


68.1±4.3 


IIZw40 


0.044±0.007 


0.044±0.006 


0.006±0.005 


0.033±0.006 


13.2±2.1 


19.0±2.8 


2.8±1.2 


21.1±3.7 


UGC4274 


0.423±0.032 


0.497±0.049 


0.105±0.012 


0.495±0.021 


23.2±1.2 


43.3±2.3 


7.7±0.9 


29.1±0.9 


IZwl8 


<0.233 






<0.116 


<0.4 






<0.1 


Mrkl450 


0.207±0.045 


0.337±0.143 


0.084±0.020 


0.112±0.017 


1.8±0.4 


2.8±0.9 


0.7±0.2 


1.6±0.2 


UM461 


<0.663 






<0.134 


<1.5 






<1.6 


Mrkl499 


0.408±0.040 


0.526±0.016 


0.156±0.011 


0.721±0.035 


5.5±0.3 


10.9±0.3 


2.3±0.2 


7.3±0.2 



Note. — Contrary to lWu et al.l j200(j|) . the wavelengths whieh we have chosen to determine the underlying continuum are fixed in 
the above measurements, which explains the minor discrepancies between these two studies. The symbol " - - ■ " indicates that no 
PAH EW measurement was possible and it is mostly due to the low SNR of the corresponding spectrum. In this case, the determination 
of the continuum is highly uncertain and could significantly affect the value of the PAH EW. For SBS0335-052E, the SNR is high 
enough, but no PAH features can be identified in its mid-IR spectrum. We do not derive upper limits for the PAH EWs at 7.7 fim and 
8.6 /im, because no reliable templates for those features are available. 

metallicity found in SBS0335-052E, but if this were true 
for the whole sample, it should also apply to the sulfur 
abundances. However, when we plot the S/H (Optical) 
against the S/H (IR) in Fig. 5b, we find that the sources 
are located on both sides of the 1:1 proportion line. Al- 
ternatively, this could be due to the difference in the 
determination of element abundance using the infrared 
and optical methods. Because only NeHI is detected 
in the optical regime, the total elemental abundance of 
neon from the optical study is heavily dependent on the 
adopted ICF. In the infrared the [NeH] line is detected 
in most of our sources and we do not use any ICF for our 
study. Thus the higher neon abundances derived from 
the infrared data compared to the optical results could 
also be due to the large uncertainty in the ICF used in 
the optical studies. Another possibility is that the tem- 
perature of the Nelll ion is lower than that of the QU I 
ion as found in some PNe (|Bernard Salas et al.l[200l ). 
Because the optical studies are based on Te(OIII) when 
calculating the ionic abundance of Nelll, that could also 
result in an underestimate of its abundance. 

Finally, we explore the variation in the ratio of Nc/S 
as derived from the infrared and the optical data, and 
plot them as a function of 0/H in Fig. 6. Interest- 
ingly, we find that with the exception of VIIZw403, in 
all BCDs, the infrared derived Ne/S ratios are higher 
than th e corresponding optical values. Ilzotov fc Thuaiil 
(11999b) indicate an average Ne/S ratio of 6.9 for the 54 
super giant HII regions they study in 50 BCDs using op- 
tical spectroscopy. The average of the optically derived 
Ne/S ratio for the nine sources in our sample is 6.5±1.8 
while the corresponding ratio from the infrared data is 
11.4±2.9, statistically higher than the optical results. As 
we mentioned earlier, no SII abundances could be de- 
termined with the Spitzer/IRS data. However, in high- 
excitation objects such as BCDs, SII does not contribute 
much ('-^10%) to the total elemental abundance, thus it 
could not account for the factor of two difference in the 
ratio of Ne/S. Furthermore, such a discrepancy is not a 
result of metallicity because we see no correlation in the 
dispersion between the infrared and optically derived ra- 
tios of Ne/S with respect to 0/H. If we divide the sources 
into a high-excitation group with sources that have a de- 
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Fig. 5. — a)The infrared derived Ne/H as a function of the 
optical Ne/H abundances. The sources shown with only the error 
bars are those where [Nell] or [SIIII] is not detected. The solid line 
represents the 1:1 line for the optical and IR derived abundances, 
b) Same as a), but for S/H (Optical) vs S/H (IR) abundances. 

tectable [OIVJA25.89 /im line and a low-excitation group 
with sources without any detection of [OIV], we still see 
no clear trend. Thus the ionization field cannot explain 
the difference in the Ne/S ratios either. Moreover, no 
correlation is found between the Ne/S ratio and the ex- 
tinction to the source. Even though the largest uncer- 
tainty in the infrared derived abundances comes from the 
uncertainty in the H/3 flux in cases where no Hu a is de- 
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Wu et al. 




Fig. 6. — The abundance ratio of Ne/S as a function of O/H. The 
infrared derived Ne/S ratios are indicated by the filled circles while 
the optical derived Ne/S ratios are shown by the open circles. The 
dashed line is the average Ne/S ratio from the infrared derived ele- 
mental abundances while the dotted line is average Ne/S ratio from 
the optical derived elemental abundan ces. The dash-dotted li ne 
represents the average Ne/S ratio from llzotov fc ThuanI {19993). 

tected, when calculating the ratio of Ne/S, the H/3 flux 
cancels out, thus our infrared determined values for Ne/S 
should be fairly reliable. One possibility for the observed 
discrepancy could be that the difference in the infrared 
and optical results is again due to the ICF adopted in 
the optical studies. 

5. METALLICITIES AND PAHS 

The observed suppression of PAH emission in_Jow 
metallicit y environm ent is still an open topic ([Maddenl 
2^00; M adden et a l. 2006; Engelbrac ht et all l200a 
Wu et al.n2006l : lO'Halloran et al.ll200QV One explana- 
tion on the absence of PAHs is that there is not enough 
carbon to form into PAH molecules in such metal-poor 
environments while an alternative possibility could be 
that the strong ionization field in low metallicity envi- 
ronments may destroy the fragile PAH rings. 

Previous studies relating the PAHs and metallicity 
were all based on the oxygen abundance measured from 
the optical, which could be associated with regions that 
are different from the ones where PAH emission origi- 
nates. Because we have obtained the infrared measured 
neon and sulfur abundances, it is particularly interest- 
ing to compare the PAH strength with our metallicity 
estimates. Several BCDs in our sample do show PAH 
emission features and an analysis of the PAH properties 
of our samp le based on IRS low resolution spectra was 
presented in lWu et"al] ()2006f ). We have obtained deeper 
spectra for some of the sources, so for consistency, we 



have re-measured all features in our sample and the re- 
sults are given in Table [5l None of the sources with 
metallicity lower than ~0.12Zq show any detectable 
PAHs. For the sources where PAHs were detected, we 
calculate the equivalent widths (EWs) of the 6.2,7.7,8.6 
and 11.2 /xm features and find that the correlation be- 
tween the PAH EWs and the metallicit y is weak. This i s 
not unexpected because as discussed in lWu et"al] (|2006D . 
the strength of PAHs in low-metallicity environments is 
a combination of creation and destruction effects. As a 
result the dependency of PAH strength on the metallicity 
alone could have a substantial scatter. 

6. CONCLUSIONS 

We have studied the neon and sulfur abundances of 
thirteen BCDs using Spitzer/IRS high- resolution spec- 
troscopy. Our analysis was based on the fine-structure 
lines and the hydrogen recombination line detected in the 
SH spectra of the IRS, combined with the radio contin- 
uum, Ha images and integrated optical spectral data in 
some cases. We find a positive correlation between the 
neon and sulfur abundances, though sulfur appears to be 
more under- abundant than neon (with respect to solar). 
The ratio of Ne/S for our sources is on average 11.4±2.9, 
which is consistent with what has been found in other 
HII regions using infrared data. However, this average 
ratio appears to be higher than the corresponding opti- 
cal value of 6.5±1.8 (in BCDs), which could be due to 
the adopted ICFs in the optical studies. When compar- 
ing the newly derived neon and sulfur abundances with 
the oxygen abundances measured from the optical lines, 
we find a good overall agreement. This indicates that 
there are few completely dust enshrouded HII regions in 
our BCDs, or if such HII regions are present, they have 
similar metallicities to the ones probed in the optical. 
Finally, the infrared derived neon and sulfur abundances 
also correlate, with some scatter, with the corresponding 
elemental abundances derived from the optical data. 
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